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ABSTRACT 
The  stress  and  strain  behavior  at  notch  tips  is  analyzed  with 
the  view  of  using  local  values  for  fatigue  calculations.   A 
perfectly  plastic,  closed  form  solution  is  used  as  a  starting 
point  and  basis  of  comparison  for  the  finite  element  analysis. 
Maximum  stress  and  strain  values  for  7075-T6  aluminum  were 
obtained  via  finite  elements  and  reflection  photoelasticity  in 
the  plastic  zone,  and  residual  stress  and  strain  values  after 
unloading  were  obtained  as  well.   A  uniaxial  model  was  used 
successfully  to  predict  most  of  the  behavior  observed. 
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A  STUDY  OF  NOTCH  FATIGUE 
PART  III:   PLASTIC  STRESS  ANALYSIS  OF  NOTCHES 

INTRODUCTION 

As  aircraft  maneuver  through  a  varied  and  complicated  load 
history,  critical  points  at  sites  of  stress  concentration  are 
locally  loaded  into  the  plastic  range.   This  is  a  good  protec- 
tion factor  for  the  structure,  for  the  yielding  provides  a  large 
load  bearing  capability  for  loads  beyond  the  limit  load.   However, 
yielding  at  the  fatigue  critical  points  also  greatly  influences 
the  fatigue  life  of  the  part,  and  calculations  of  fatigue  must 
include  these  effects  in  order  to  be  accurate. 
RESIDUAL  STRESS 

After  a  load  cycle  is  applied  which  takes  the  ligaments  at 
the  notch  tip  into  the  yield  zone,  those  ligaments  have  a  perma- 
nent set,  which  makes  them  longer  than  their  original  length. 
During  the  unloading  portion  of  the  cycle,  the  surrounding  un- 
yielded  material  attempts  to  push  the  elongated  ligaments  back 
to  their  original  length  and  a  position  of  equilibrium  is  attained 
with  the  ligaments  in  compression;  thus,  a  compressive  residual 
stress  is  produced  in  the  ligaments.   The  ligaments  are  still 
longer  than  they  were  in  the  unyielded  state  and  the  notch  tip 
region  possesses  a  zone  of  compressive  residual  stress  but  posi- 
tive residual  strain.   Both  of  these  must  be  accounted  for  in 
a  cycle  by  cycle  summation  of  damage  for  fatigue  life  monitoring. 
STRESS  CONCENTRATION  FACTOR 

Tensile  yielding  at  the  notch  tip  causes  a  reduction  in 
the  stress  concentration  factor.   Thinking  of  the  ligament  at 
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the  notch  tip  as  a  miniature  uniaxial  tensile  specimen,  as 
the  stress-strain  curve  for  the  ligament  bends  over  with  yield- 
ing, the  ratio  of  far-field  stress  to  the  local  plastic  stress 
produced  is  larger  than  for  the  elastic  condition  and  this  rela- 
tionship reduces  the  SCF.   This  change  in  SCF  can  have  signi- 
ficant effect  on  fatigue  life  predictions  for  those  cycles 
applied  after  a  cycle  which  causes  yielding. 
PHOTOELASTIC  RAMIFICATIONS 

After  yielding  occurs  in  the  specimen  material,  stresses 
and  strains  in  the  coating  material  undergo  a  shift  relative  to 
one  another  because  the  coating  strains  are  constrained  to  be 
the  same  as  the  yielded  metal  so  that  the  classical  relation- 
ships for  stress  in  terms  of  fringe  number  no  longer  hold.   In 
this  report  strains  are  obtained  photoelastically ,  and  the 
stresses  are  calculated  using  plastic  stress-strain  laws.   The 
coating,  however,  acts  as  a  permanent  repository  of  the  residual 
strain  for  whatever  loading  may  follow.   To  obtain  the  stresses, 
a  model  must  be  formulated,  which  gives  the  new  stress-strain 
relation,  where  the  zero  stress  point  doesn't  coincide  any  longer 
with  the  zero  strain  point  but  the  two  have  been  shifted  apart. 

All  three  of  these  phenomena  will  be  explored  in  detail  in 
this  report  with  supporting  experimental  data  to  quantify  the 
behavior. 

PLANE  STRESS  SLIPLINE  THEORY 

To  begin  exploration  into  the  nonlinear  domain,  some  ana- 
lytical methods  were  first  searched  for,  and  the  only  candidate 
seemed  to  be  slipline  theory.   The  material  model  is  highly 


viealized,  but  with  photoelasticity  and  FEA,  it  is  possible 

to  assess  its  usefulness.   It  is  an  elegant  theory,  and  it 

proved  very  useful  in  checking  out  the  nonlinear  portion  of  the 

ADINA  Program.   The  development  here  follows  that  given  by 

Kachnov 

GOVERNING  EQUATIONS 

In  the  plastic  region,  slipline  theory  assumes  that  the 
stress  state  is  on  the  original  yield  surface,  or  in  other  words, 
that  the  material  is  perfectly  plastic.   Using  the  Von  Mises 
yield  criterion,  this  requires  that 


(Tj'-cr.cr,  +oV  =  3lc* 


(i: 


where   k  =  yield  limit  in  pure  shear.   The  uniaxial  yield  stress 
is  Gy  :  i5  ^    .   A  parametric  form  of  this  equation  was 
developed  by  Sokolovskii, 


GTr  2kdos  (to-ir)      CTt=  Z\cC*>{»>  +  TL) 


2) 


which  satisfies  (1)  identically  for  all  values  of  u  . 

The  parametric  function  u)(x,y)   is  related  to  the  mean  pressure, 

which  is  given  by 


(T-    *(<r>  +  ^  +  <r*) 


5V--       -       -  (3 

Substituting  Equation  (2)  into  Equation  (3)  and  setting  0^  =  ° 
for  plane  stress,  we  find 

Co*  VJ     -  k'  (4) 


To  obtain  the  expression  for  equilibrium  equations  in  terms 
of   u)  ,  the  transformation  equations  are  first  used 

£  =  i  (<j7 + d)  *  £  (j;  -crt)  c*  2/ 


(5a) 


G"y  =  i  (<r,  ♦  <rt)  -  ^(<r,  -<£)  Co  2^  (5b) 

where   $   is  the  angle  between  the  first  principal  direction  and 

the  x   axis. 

Substituting  Equations  (2)  into  (5), 


(7-  uf^^0*60  -   ScrVlO  Cos  *<*] 


(6a) 

(6b) 
(6c) 


Substituting  the  stress  components  into  the  differential  equations 
of  equilibrium  simplified  to  plane  stress,  we  obtain  a  system  of 
two  equations  for  the  two  unknown  functions   u)(x,y)   and   4>(x,y)  . 
After  some  rearrangement  to  simplify  later  operations,  the  two 
equations  become 

ftf5StAio  Cos  Z^- Co*  to)  ^lO   4  ^SUco  ^i*l4c>*> 

-Z.SUuJ  H  -  o 

(ftS\KU>Sl**/W    -(fF5Uu0C0S2^H-C0$lo)^0  (7b) 


METHOD  OF  CHARACTERISTICS 

Seeking  a  solution  by  the  method  of  characteristics,  we 
suppose  that  along  some  line  in  the  structure  given  by   x  =  x(s, 
and  y  =  y(s)  ,  the  functions  u>   =  oo(s)   and   4>  =  <$>{s)      are 
given.   Then  along  that  line 


2s      3x7  T?>      ^y  35 


(8a) 


du>  .  2">  dx  4.  ^  <^y 

Zi  "  ~?x*   Js      "?y   3"s  (8b) 

Equations  (7)  and  (8)  possess  characteristics  whose  defining 
differential  equations  are  found  by  setting  the  determinant  of 
the  coefficients  to  zero.   Reducing  the  determinant  we  obtain 


"3k"    V?  su*oeoii/-cei 


lO 


(9) 


In  order  for  the  characteristics  to  be  real,  the  discri- 
minant must  be  positive  semi  definite, 

3-4Co»*«j2o  (10) 

or  cO  *   1L 

Therefore,  because  of  the  nature  of  the  cosine  function,  the 
region  of  hyperbolic  characteristics  lies  between  the  limits 
ft"  ^  UO  <m   Sjtt"  .   From  Equation  (10) 

Z~     '  u     ,      <  * 

(2o<>  oo  <  —  (11) 

and  from  Equation  (4)  we  see  that  this  physically  corresponds  to 

CT5  k.  (12) 


which  means  that  the  mean  pressure  must  be  less  than  the 

maximum  shear  stress  in  order  to  have  a  region  where  the 

equations  are  hyperbolic  and  characteristics  exist. 

In  the  aoolication  of  Kramer's  Rule  to  solve  the  set  of 

simultaneous  ecmations  for  dx,  a  numerator  determinant  is  formed 

ds 
Setting  this  determinant  to  zero  gives  the  condition  along  one 

set  of  characteristics,  which  reduces  to 


ft   +  d    =   CoA^taJvt 


(13) 


where 


(14) 


Along  the  other  set  of  characteristics,  obtained  by  setting 

the  determinant  for  ■%*-   =  o  ,  we  find 

ds 


jx  -4  -  Co**t<u*.t 


(15) 


CIRCULAR  BOUNDARIES 

Consider  a  portion  of  a  circular  boundary  as  shown  in 
Figure  1. 


/3  characteristic 


o<.  characteristic 


FIGURE  1   SCHEMATIC  OF  CHARACTERISTICS 


Stress  conditions  at  point  P  will  only  be  influenced  by  condi- 
tions along  the  portion  of  the  boundary  between  A  and  B,  and 
if  the  boundary  changes  contour  outside  of  AB  it  will  not  affect 
the  stresses  at  P.   More  specifically,  the  stresses  at  P  for  a 
notch  and  for  an  axi symmetric  hole  would  be  the  same,  so  long 
as  all  of  the  material  between  A  and  B  were  yielded  and  the 
boundary  is  stress  free.   At   r  =  a  all  along  the  yielded 
boundary,  x .  . Q  =  0   and   a   =  0  .   By  the  Von  Mises  yield  con- 
dition in  Equation  (1)  ,  (jL  z.   |  3  k-   ;  the  mean  pressure  is 
Gr  —-  ;  and  the  maximum  shear  stress  is   *«%•.*  =  1—  K.  ;  thus, 
the  condition  is  satisfied  for  characteristics.   From  Equation 
(4)   oj  =  -^   at  the  boundary   r  =  a  . 

From  Equation  (2)  the  normal  stresses  in  the  yielded  region 
are  given  by 


(T    =  ZJL   Co*  (w  '  £) 


(16a) 
(16b) 


STRESS  DISTRIBUTIONS 

To  relate   w  to   r  ,    the  axisymmetric  form  of  the  equilibrium 
equations  is  used. 

-r-  + "  ^  (17) 

<:\r      r 

Substituting  Equations  (16)  into  (17)  and  simplifying,  a  relation 
between   u   and   r   is  obtained. 


({*  *  Cot  l»j)A*>  +  *fk*  ^o 


(18 


Integrating, 


(19) 


where  C   is  a  constant  of  integration. 

As  shown  before,  on  the  boundary  at  r  =  a  ,  co  =  ■=-  .   Using 
this  point  to  evaluate   C  . 


z 


20) 


r  - 


£  a. 

£    5U^  (21) 


TT 


which 


Characteristics  will  extend  out  into  the  body  to  a>  =  ^  , 
corresponds  to   r  =  2.07a  .   Within  that  region  this  solution 
will  apply  to  the  notch. 

Equation  (21)  has  been  evaluated  for  the  notch  in  Figure  2 
given  by   r  =  1.8125  in.  and  d  =  2.875  in.  out  to  a  distance 
of   r  =  3.75  in.,  which  corresponds  to  the  limit  of  applicability 
of  the  solution.   The  results  are  listed  in  Table  I. 

The  domain  of  applicability  of  the  solution  is  defined  by 
the  free  boundary  at  the  edge  of  the  notch  and  the  two  char- 
acteristics given  by  Equations  (13)  and  (15).   Since  the  solu- 
tion for  the  notch  coincides  with  the  axisymmetric  solution  in 
this  range,  a        and   a  Q   are  principal  stresses  and   <j>  =  9  . 
Thus, 

(22) 


O  i  &  -    Co**t<v-^t 


LARGE    SPECIMEN    GEOMETRIES 


96" 


p 

d 

23.5' 


Nominal    Ky  2.6  3.8 

r  (in )        1,8125       0.625 
d  (in  )        2.875         3.938 


Reduced 
Cross-section 
(in2) 


l.42in2       1.25in2 


7075-T6    Aluminum 
Thickness     0.080  in 


FIGURE  ?.  NOTCH  GEOMETRY 


TABLE  I 

NOTCH  STRESS  DISTRIBUTION  BY  SLIPLINE  THEORY 

r         r/a  co  a    //Ik  a  .//3k 

1.8125       1.0000  tt/3     =  1.047  0.0000  1.0000 

2.1645       1.1942  tt/3. 49  =   .900  0.1689  1.0740 

2.4698       1.3626  tt/3. 93  =   .799  0.2832  1.1113 

2.8381       1.5658  it/4.49     .700  0.3925  1.1369 

3.3057       1.8238  ir/5.24     .600  0.4989  1.1514 

3.7506       2.0690  tt/6        .524  0.5774  1.1547 
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At  the  outer  extremity  of  the  zone,  where   ui  =  tt/6   and 
r  =  3.7506  ,  the  characteristics  intersect  and  coincide.   From 
the  definition  of   8   given  in  Equation  (14),  it  is  zero  at 
u)  =  tt/6   for  all   6   .   Concerning  ourselves  with  the  point  P 
in  Figure  1,  which  corresponds  to  the  symmetry  line  of  the 
notch  where   9=0,  Equation  (22)  yields  the  value  of  the 
constant  as  zero.   We  therefore  are  looking  for  the  character- 
istics defined  by 

Jl  t  G     -    &  (23) 

Evaluating   ft   from  Equation  (14)  by  Simpson's  rule  at  all  of 
the    r   values  in  Table  I,  and  obtaining  the  corresponding   9 
coordinates  from  Equation  (23),  the  trajectory  of  the  character- 
istics can  be  found.   They  have  been  tabulated  in  Table  II. 

A  plot  of  these  contours  has  been  made  on  the  tip  of  the 
notch  in  Figure  3.   This  is  not  a  drawing  of  the  plastic  zone, 
but  is  the  boundary  of  applicability  of  the  slipline  solution. 

This  slipline  method  provides  about  the  only  analytical 
means  of  handling  the  nonlinear,  plastic  flow  at  a  notch  tip. 
Although  it  is  limited  to  rigid,  perfectly  plastic  material 
models  and  the  extent  of  yielding  is  unknown;  it  does  constitute 
a  valid  check  on  the  FEA  and  later  these  results  will  be  compared 
to  those  attained  for  more  realistic  material  models. 
COMPARISON  WITH  FEA 

The  finite  element  grid  used  for  the  elastic  analysis  re- 

(2) 
ported  in  Part  II     was  employed  with  ADINA's  bilinear  material 

model  such  that  it  represented  a  rigid  perfectly  plastic  mate- 

rial.   A  Young's  Modulus  of  10    psi  was  used  to  model  perfect 

11 


TABLE  II 


COORDINATES  OF  BOUNDING  CHARACTERISTICS 


w 


1.8125 

1.047 

+ 

24.1° 

2.1645 

0.900 

+ 

17.3° 

2.4698 

0.799 

+ 

12.7° 

2.8381 

0.700 

+ 

8.1° 

3.3057 

0.600 

+ 

3.5° 

3.7506 

0.524 

0° 
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rigidity,  and  the  strain  hardening  modulus,  E  ,  was  set  equal 
to  zero.  Poisson's  ratio  was  entered  as  0.4999999,  which  is  as 
close  to  incompressibility  as  the  program  would  allow. 

Figure  4  illustrates  the  FEA  results  obtained  and  compares 
them  to  the  slipline  solution.   The  stress  values  are  normalized 
to  the  yield  stress.   Very  good  correlation  was  obtained  for  the 

a .,   distribution  but  the   a    values  vary  somewhat.   This  seems 

9  r  J 

to  be  par  for  this  problem  in  that  it  seems  no  matter  what  method 
is  used  to  find  the  stresses,  the  greater  variation  is  always  in 
a   .   The  comparisons  are  viewed,  however,  as  substantiations  of 
one  another  and  constitute  a  baseline  to  which  other  analyses 
can  be  referenced. 

The  growth  of  the  plastic  zone,  which  is  not  obtainable  in 
the  slipline   theory,  is  given  by  FEA  in  Figures  5  to  7 . 

We  now  explore  a  more  realistic  material  model  with  FEA  and 
with  photoelasticity .   In  order  to  do  that,  it  is  necessary  to 
develop  a  nonlinear  stress-strain  law. 

PLASTIC  STRESS-STRAIN  RELATIONSHIP 
When  the  proportional  limit  is  exceeded  at  the  notch  tip,  . 
stresses  can  no  longer  be  computed  directly  from  the  photoelastic 
fringe  data.   A  plastic  stress-strain  relationship  must  be  used  to 
calculate  stresses  from  measured  strain  values.   From  Druckers's 
Postulate     it  has  been  shown  that  the  plastic  strain  increment  is 


normal  to  the  yield  function  in  stress  ;  i.e. 

(24) 


14 


ui      .. 


in 
u 


5fr     in 

CE 


Q. 

I 


U 
K 

ui 

Lu 

ac 

CL 
I 


n* 


03 


si 


si 


si 


D/D '  UUUHJ  553H1S 


Si 

x 


l/i 
ri 


ui 


fM 


z 
m 


u 

M 


CL 


Ln 


a 

si 


15 


FIGURE  5 


RIGID-PERFECTLY-PLASTIC   INITIAL  PLASTIC  ZONE 
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FIGURE  6 
RIGID-PERFECTLY-PLASTIC  INTERMEDIATE  PLASTIC  ZONE 
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FIGURE  7 


RIGID-PERFECTLY-PLASTIC   FINAL  PLASTIC  ZONE 
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where   deij^  =  Plastic  strain  increment 

f(aij)=  0  describes  the  yield  function 

dN=  An  undetermined  function  of  stress,  strain  and 
history  of  loading. 
The  Von  Mises,  or  distortional  energy,  yield  criterion  is 
good  for  aluminum.   In  terms  of  principal  stress  invariants, 
yielding  occurs  when   J~   reaches  a  given  value. 


Jt  =  t\«rt-<rS+&  -^)l*  to-o^y  }  -_  (Vst~t 

where   J~  =  Second  deviatoric  stress  invariant 
a .    =  Principal  Stresses 

The   f   is  defined  as 

Differentiating  as  indicated  in  Equation  (24). 


(25) 


(26) 


27a) 


(27b) 


(27c) 


* 


To  determine   dA  ,  an  equivalent  uniaxial  stress,  a      ,  is 
defined  to  represent  the  behavior  of  a  triaxial  stress  state 
with  a  single  component. 
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(J*  5  F(^;crl(<r5) 


(28 


such  that  in  uniaxial  tension, 


r  (61.0,0)  =  cr; 


(29 


and  at  yield  when   f  =  0 


G*--C, 


(30) 


when   f  «  0  /  J-  -  Const  ,  so  it  is  convenient  to  let 


ct*=f(jO 


The  simplest  function  is  customarily  selected 


(31) 


32 


where  the  square  root  is  chosen  to  make  the  expression  dimension- 
allv  correct.   For  uniaxial  tension, 


F(<r,,o,o)*.M2(rl'-}t 


(33) 


To  satisfy  the  condition  of  Equation  (29) , 


A  =  /3 


and  the  equivalent  stress  becomes 


f.  J=  f(<r,-<ro*+  fa-«tf*  C<r'-^)t]t 


34) 
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A  generalized  plastic  strain  increment  function  is  defined 
such  that  the  product  of  the  two  is  work. 


After  considerable  manipulation, 


35 


36 


Inserting  Equations  (27)  into  Equation  (36)  and  simplifying  with 
Equation  (34) ,  it  can  be  shown  that 


(37 


So  that  the  stress-strain  law  now  becomes 

(T*  (38a) 


38b) 


<n;-[c*-i  &+«*)]  ip 


clfc* 

+  (38c) 


*  P  - 


Simplifying  Equation  (38a)  to  uniaxial  tension,   de,^  =  de 

*  .        .        d     * 
because  a      =  a,  ;  thus,  upon  integration  ,  e,1-  =  e   ,  since 

*  P        P        P 

e   =0   when  e,       =  e2   =  £3   =  °  •   This  leads  to  a  relationship 

*  * 
between   a    and   e 

It  is  postulated  that  if  in  uniaxial  tension,  for  loading 

t[-      F(0  (39) 
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then  a  reasonable  generalization  is 


£*-  f(^) 


40) 


For  instance,  using  the  Ramburg-Osgood  expression  for  uniaxial 
tensile  behavior, 


*'.-/*(%) 


(41) 


then  it  is  proposed  that 


•/*  ar 


(42) 


and 


*.  W-l 


at*c  *£(§]    <kr* 


E^e 


43) 


The  stress-strain  expression  now  has  the  form 


*t!:k-t{w)Uu<£) 


%\  V\-2- 


At 


(44) 


For  the  notch  specimens  the  loading  was  proportional;  i.e. 
Under  these  conditions  the  plastic  strain  becomes 


45) 


m;-.  w^KIt)  cW* 


(46) 


where   y  =  Const  =  1  -  1/2 (a,  +  a-) 
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Also  generalized  stress  from  Equation  (34)  becomes, 


(47) 


where   5  =  Const  = 


-if  (| -«i,»)  +  (cl%-dtf  ■*(»-<)]*' 


Equation  (46)  then  simplifies  to 


48) 


This  is  now  integrable,  and   e,p   can  be  found  explicitely 


<-  *~V  (!}" 


(49 


The  integration  constant  is  zero  from  Equation  (41) . 

Eliminating  y      and   5   using  Equations  (46)  and  (47). -the  plastic 

strain  can  be  expressed  in  terms  of  stresses. 


'=[S-i(O<r0]4(f) 


*,'  = 


*\  A-\ 


(50) 


The  elastic  and  plastic  parts  of  the  strain  are  assumed  additive 


e,  =  c  +  t; 


(51) 


and  the  total  strain  becomes 


e»s|[6;-0ffi*«*)]+ 4  (£)*  '[or-^M) 


(52) 
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For  plane  stress,  a      =  0  ,  and  both  of  the  other  components  can 
be  obtained  by  permutation. 


(53a) 


(53b) 


53c' 


This  stress-strain  law  will  be  used  to  find  stresses  from 
photoelastic  strain  data,  and  all  other  analysis  of  the  plastic 
zcne  around  the  notch  except  for  the  work  done  with  the  rigid  - 
perfectly  plastic  model,  which  was  used  for  the  slip-line  theory. 
Since  the  strains  in  the  coating  are  constrained  to  be  the  same 
as  the  specimen,  the  equations  relating  strains  and  fringe  will 
remain  valid  into  the  plastic  domain.   With   e,   and   e  „   deter- 
mined from  the  photoelastic  data,  a  program  for  the  solution  of 
simultaneous,  nonlinear,  algebraic  equations  was  used  to  solve 
Equations  (53a)  and  (53b)  for   a,   and   a_  .   With   a,   and  a 2  , 
e,   could  be  found  from  Equation  (53c)  if  desired. 
PLASTIC  MEASUREMENTS  AT  THE  NOTCH 

The  shallow  notch  specimen  shown  in  Figure  2  was  fabricated 
from  7075-T6  and  loaded  into  the  plastic  range.   Normal  and 
oblique  incidence  fringes  were  recorded  in  1/16"  increments  on 
the  line  of  symmetry  extending  from  the  notch  tip  beginning  at 
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1/8  inches  and  extending  to  1/2  inches.   Five  different  readings 
were  taken  at  each  location  while  the  specimen  was  under  70,000 
lbs.  of  load.   The  readings,  along  with  the  mean  and  the  standard 
deviation,  are  recorded  in  Table  III  and  Table  IV. 

With  five  samples  taken,  data  falling  outside  of  an  expected 
value  of  l/2n  ,  or  1/10,  were  rejected.   For  a  normal  distribu- 
tion this  is   4-1. 645a   from  the  mean.   The  data  were  examined 
for  this  tolerance  and  one  oblique  reading  at   x  =  .1875  in.  fell 
outside  the  limits  and  was  discarded.   A  new  average  was  computed 
at  this  location  of  7.71. 
EXTRAPOLATION 

A  linear  regression  ignoring  the  sign  was  run  for  both  sets 
of  data  using  log  (log  N)  vs.  x  as  the  function  that  is  nearly 
linear.   Results  of  the  regression  for  normal  incidence  gave 

log  (log  N)  =  .03627  -  .  2594x  (54) 

with  R  =  .9997.  At  the  edge  of  the  hole,  x  =  0  ,  N  =  -  12 . 22 . 
(The  negative  sign  was  again  affixed  after  the  value  was  found. 
For  oblique  incidence, 

log  (log  N  )  =  .01588  -  .  3566x  (55) 

with   R  =  .99  84.   At  the  edge  of  the  hole  ,  NQ  =  -  10.89  5. 
POISSON'S  RATIO  iMISMATCH 

This  extrapolation  was  made  from  a  region  of  the  sheet 
where  the  Poisson's  Ratio  of  the  metal  is  0.5  due  to  yielding. 
From  Lindsey^4)  the  effective  Poisson's  Ratio  in  the  coating 
material  is  given  by 
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TABLE  III   FRINGE  VALUES 
SHALLOW  NOTCH  TEST  -  NORMAL  INCIDENCE  -  PLASTIC  LOAD 


AUG       S.D. 


,5000  -  6.36  -  6.39  -  6.40  -  6.42  -  6.39  -  6.392  .0217 

,4375  -  6.83  -  6.83  -  6.85  -  6.87  -  6.84  -  6.844  .0167 

,3750  -  7.44  -  7.44  -  7.41  -  7.41  -  7.42  -  7.424  .0152 

3125  -  7.98  -  7.99  -  8.03  -  8.08  -  8.06  -  8.028  .0432 

2500  -  8.65  -  8.64  -  8.65  -  8.67  -  8.67  -  8.656  .0134 

1875  -  9.29  -  9.33  -  9.38  -  9.37  -  9.35  -  9.344  .0358 

1250  -  10.11  -  10.20  -  10.19  -  10.16  -  10.23  -  10.177  .0455 


TABLE    IV      FRINGE    VALUES 
SHALLOW    NOTCH    TEST    -    OBLIQUE    INCIDENCE    -    PLASTIC    LOAD 


AUG  S.D. 


500  -  4.79  -  4.79  -  4.79  -  4.79  -  4.79  -  4.790  0 

4375  -  5.41  -  5.39  -  5.36  -  5.32  -  5.32  -  5.360  .0406 

3750  -  5.85  -  5.85  -  5.84  -  5.84  -  5.86  -  5.848  .0084 

,3125  -  6.41  -  6.40  -  6.41  -  6.40  -  6.37  -  6.398  .0164 

2500  -  7.06  -  6.97  -  7.03  -  6.89  -  6.97  -  6.984  .0654 

1875  -  7.79  -  7.70  -  7.71  -  7.73  -  7.70  -  7.726  .0378 

1250  -  0.62  -  8.57  -  8.55  -  8.58  -  8.58  -  8.58  .0255 


26 


4797  -.1207 


t 


-o- 


Cosh  C  Sir) 


■] 


.992.3  "  -C4S5  Costrv  jzE±} 
LCosk    (Str) 


56) 


Where  x*  is  measured  from  an  interior  point  far  from  the  edge  of 

the  coating.   To  relate  it  to  x  which  is  measured  inward  from  the 

notch  edge,  let  x*  =  10  ,   -  x.   For  photoelastic  measurements  at 

points  where  x>  .1250  in,  Poisson's  Ratio  is  essentially  constant 

The  coating  being  used  was  he  =  .0382  ;thus  x  =  .1250  corresponds 

to  x  -  3.2723h   and  x*  =  6.7277h  .   Substituting  this  value  into 
c  c 

(56)  gives  v  =  .4828. 

However  when  extrapolating  back  to  x  =  0  ,  the  extrapolation  is 

from  data  for  which  the  value  of    is  .4828  but  the  actual  value 

(4) 
there  from  (56)  is  .3793.   Lindsey    worked  out  the  correction 

for  the  extraDolation : 


H, •(«["],.. 


V=.3795 


8ZS 


(57) 


The  same  correction  applies  to  the  oblique  incidence  fringes; 
therefore,  the  fringe  value  at  the  edge  of  the  notch  is 

N:  (.9Soz)(-l2.«}=  -11-37     Ne--(.93oz)(-.o.9o)---lo.l3  (58 

REINFORCEMENT 

A  slight  reinforcement  of  the  thin  metal  sheet  by  the  coat- 
ing material  is  present,  and  this  small  correction  factor  was 
develoDed  in  Part  I. 


C-  C 


r        htEcO+0,) 


h4Es(l+Ot) 
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(£,*-£«) 


(59) 


In  the  elastic  domain,  the  factor  was  1.0162,  but  in  the  plastic 
domain  it  is  1.0182.   This  correction  only  need  be  made  if  com- 
parisons are  being  made  with  calculations  or  measurements  being 
made  on  an  uncoated  part,  and  this  will  be  the  case  when 
comparisons  are  made  later  with  finite  element  solutions.   The 
finally  corrected  values  of  the  fringes  are  given  in  Table  V. 
STRAINS 

The  strain  expressions  reflecting  the  calibrated  parameters 

of  the  polariscope  in  the  Aeronautical  Engineering  laboratory  at 
(4) 


NPS  are 


(60a! 


60b! 


c*     o+a)L  J 

where  z      -  e   =  F  N  and  e   may  or  may  not  be  the  larger 
x    y    e       x   J  J  r 

principal  strain.   At  the  edge  of  the  hole  where  v   =  .3793 

sK=.  -l.5798*.io'*N  +  2.429i*io'4N« 

At  the  next  station  and  all  remaining  stations,  x£0.1250in., 
Poisson's  Ratio  is  0.4828  and  the  strain  expressions  are  given  by 

U*  (62a) 


(61a) 


g*-  -2.79oi*io-*M  +  I.8&MXIO  Me 


(62b) 
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TABLE  V   CORRECTED  FRINGE  VALUES 


x  N  N, 


0.0000 

- 

11.58 

-  10.31 

.1250 

- 

10.36 

-   8.74 

.1875 

- 

9.51 

-   7.85 

.2500 

- 

8.81 

-   7.11 

.3125 

- 

8.17 

-   6.51 

.3750 

- 

7.56 

-   5.95 

.4375 

- 

6.97 

-   5.46 

.5000 

_ 

6.51 

-   4.88 
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Using  the  fringe  values  of  Table  V  in  Equations  (61)  and  (62) , 

the  strain  distribution  is  obtained  in  Table  VI. 

STRESSES 

The  nonlinear  constitutive  relations  developed  in  Equations 
(53a)  and  53b)  were  used  to  find  stresses.   Since  the  equations 
cannot  be  solved  explicitly  for  stress,  a  computer  program  was 
used  to  solve  the  nonlinear  set  of  equations  simultaneously. 
Results  in  the  form  of  a  stress  distribution  are  given  in 

Table  VII,  where  the  following  values  of  the  parameters  were 

(2) 
used  as  given  by  Lindsey    in  Part  II. 

E  =  10.12  x  106  psi 
v  =  0.5 

6  =  1.479  x  1043 
n  =  21.58 

The  a        value  at  x  =  0  should  be  zero  of  course;  however 
it  is  less  than  1%  of  the  other  stress  component  and  represents 
a  small  experimental  error.   To  continue  the  investigation  of 
the  plastic  region  at  the  notch,  a  finite  element  analysis  was 
conducted  to  be  used  in  conjunction  with  the  photoelastic  results 

FINITE  ELEMENT  ANALYSIS 

The  ADINA  Program  was  used  as  described  in  Part  II.   The 
grids  for  the  notch  are  shown  in  Figures  8  and  9.   The  analysis 
for  loading  into  the  plastic  region  of  the  7075-T6  aluminum  was 
made  using  a  bilinear  material  model  mentioned  earlier.   The 
values  required  by  ADINA  were  determined  graphically  from  experi- 
mental data  shown  in  Figure  10  and  the  bilinear  fit  is  also 
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TABLE  VI   PLASTIC  STRAIN  DISTRIBUTION 


3  2 

x  e   x  10  ex  10 

x  v 


0.0  -  6.7551  1.3672 

0.1250  -  5.8288  1.2446 

0.1875  -  5.0249  1.1751 

0.2500  -  4.3496  1.1191 

0.3125  -  3.8764  1.0535 

0.3750  -  3.4477  .9888 

0.4375  -  3.1304  .9165 

0.5000  -  2.5101  .8974 
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TABLE  VII   STRESS  DISTRIBUTION 


X 

ax  (psi) 

°y  (psi) 

0 

629 

79,945 

.1250 

3,327 

80,441 

.1875 

7,533 

81,813 

.2500 

11,466 

82,975 

.3125 

13,447 

83,102 

.3750 

15,214 

83,012 

.4375 

15,696 

82,145 

.5000 

21,474 

83,847 
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FIGURE   8 
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FIGURE  9    FINE  MESH  FOR  SHALLOW  NOTCH 
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FIGURE  10 
707S-TE  HLUMINUM   5TRE55-5TRH !N  CURVE 
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shown.   The  loads  used  in  the  program  were  selected  to  match 
the  experiments. 

The  results  of  FEA  for  a  local  yield  producing  load  of 
70,000  lbs  are  presented  along  with  the  experimental  results  in 
Figure  11.   The  a        results  compare  well/  matching  in  both 
magnitude  and  trend  with  the  experimental  data.   In  all  cases 
the  FEA  determined  the  peak   afl   stress  to  occur  near  the  yield 
boundary,  and  the  gradient  of  the   a    stress  to  fall  off  dra- 
matically in  the  plastic  zone.   This  characteristic  behavior  of 
the   a    stress  was  reported  by  other  investigators  [References 
5  and  6]  using  FEA  on  2024-T3  aluminum.   Plane  elastic-plastic 
stress  distributions  reported  by  Frocht  [Reference  7]  show 
similar  trends.   The  experimental  data  also  shows  a  marked  change 
in  the  gradient  of   a.   stress  within  the  plastic  region.   The 
growth  of  this  plastic  region  is  approximated  by  using  the  FEA 
results  for  this  notch  loaded  at  60,000,  65,000  and  70,000  lbs. 
These  results  are  shown  in  Figures  12  through  14.   Experimental 

data  for  the   aQ   stress  distribution  matches  the  FEA  results 

y 

reasonably  close. 

The  compatibility  of  the  data  is  satisfying  because  the 
finite  element  mesh  cannot  be  made  too  fine  without  becoming 
unwieldly,  while  the  photoelastic  data  is  continuous.   Also 

because  of  the  averaging  techniques  used  to  extrapolate  the 
finite  element  results,  the  stress  values  in  regions  of  high 
gradients  are  suspect.   On  the  other  hand,  the  finite  element 
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FIGURE    II 
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FIGURE  12 


SHALLOW  NOTCH   60,000  LB  LOAD  PLASTIC  ZONE 
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FIGURE  13 


SHALLOW  NOTCH   65,000  LB  LOAD  PLASTIC  ZONE 
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FIGURE  11+ 


SHALLOW  NOTCH   70,000  LB  LOAD  PLASTIC  ZONE 
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solution  is  readily  obtained  over  the  whole  elastic-plastic 
region  and  that  would  be  very  laborious  using  photoelascity 
by  running  them  both,  reliable  results  seem  to  be  reasonably 
certain. 
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STRESS  AND  STRAIN  CONCENTRATION  FACTORS 

Turning  attention  to  the  notch  boundary  only,  a  develop- 
ment is  made  whereby  the  local  hoop  stress  can  be  calculated. 

In  the  elastic  range,  stress  and  hoop  stress  can  be  calculated 
values  to  nominal  values  are  constant  and  equal. 

K   =  K   =  K_ 

a     e     T 

In  this  definition  the  stress  concentration  factor  was  defined 


as 

local  stress  over  smallest  cross-section  stress.   In  the  yielded 
domain,  these  elastic  relationships  chancre,   K   and   K   are 

no  longer  equal  and  their  values  depend  upon  the  amount  of 

plastic  yielding  that  takes  place  at  the  notch-tip. 

LOADING  INTO  THE  PLASTIC  RANGE 

To  show  these  relationships,  Figure  15a  has  been  constructed 
It  is  for  a  most  general  case  showing  local  conditions  versus 
far-field  conditions  of  both  stress  and  strain.   The  stress- 
strain  curve  which  corresponds  to  the  loading  cycle  in  Figure 
15a  is  shown  in  Figure  15b  with  all  corresponding  points  labeled 
the  same. 

Loading  begins  at  1  with  yielding  occurring  at  2;  so  from 
1  to  2, 

f--^s--^Ts  (S3) 

Loading  continues  to  point  4  where  unloading  begins  and  the 

cycle  is  completed  at  point  5,  where  there  is  positive  residual 

strain  and  negative  residual  stress. 

42 


E«4 


E*3=  °i 


Fig.  15a   Stress  -  Stress 
Diagram 


Fiq.  15b   Stress  -  Strain  Diagrams 
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UNLOADING 

During  unloading,  and  for  subsequent  loading  cycles  up  to 
yield,  the  stress  concentration  factor  is  elastic  and  constant, 
but  different  from  the  original  Km.   This  is  shown  in  Figure  15a 

_  6~4  -  (T«. 


KT  = 


S4  (65) 

Since  Xf  -  G+ / <,A  , 

— T-  (66) 

During  the  same  portion  of  the  cycle  the  unloading  strain  concen- 
tration factor  from  Figure  15a  is 


^4  (67) 

Since  \C<  -   ^4/e+  ,  and   S4  -  t-^4-  j 


\Cr     =    Kf 


BZ.1L 


S4  (68) 

From  Figure    15b, 


69) 


Solving  for  a    and  inserting  it  into  Equation  (65) 


s4      s4 


(70) 
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Simplifying  by  summing  the  first  two  terms  to  zero  and  using 

4  =  e  4, 


the  elastic  relation  S A    =  < 


S4.  (71) 


Comparing  equation  (71)  with  equation  (68) 

(72) 
This  result  for  the  uniaxial  model  has  some  rather  far-reaching 
implications.   One  of  them  is  that  the  strain  concentration  factor 
is  the  same  as  the  elastic  stress  concentration  factor  even  into 
the  plastic  domain.   This  is  shown  by  returning  to  equation  (71) . 
From  Figure  15a  for  parallel  unloading  curves  of  stress  and  strain, 


EL* -C*   =  EtA  ~<U 


Solving  for  Ee:_.  and  inserting  it  into  (71) 

-  V 


(73) 


Kt  '-   Kfc  -(0\-CT4-H£d4) 


(74) 


Since  Ee,  =  K„  and  a4  "  aR  =  KT ' >    equation  (74)  reduces  to 

Kt  -   ^T  (75) 

which  establishes  Ke  in  the  plastic  region  to  be  a  constant. 
This  is  a  very  important  result  which  can  be  checked  experimentally 
Inserting  expressions  of  stress  and  strain  into  (75) 

C4       5 4      S4 
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Since  a~    is  an  elastic  value  and  related  to  L  through  Young's 
modulus,  we  may  conclude 

(77) 
This  simplification,  with  the  others  already  noted  reduces 
Figure  15  to  that  appearing  in  Figure  16.   In  the  second  appli- 
cation of  load,  K  «  K'  and  K   =  K'  until  the  previous 

El  G  1 

highest  stress  and  strain  values  are  exceeded. 
RELATING  LOCAL  STRESS  TO  FAR-FIELD  STRESS 

After  vielding  has  taken  place,  we  desire  to  relate  K   to 

a 

K   so  that  local  stresses  can  be  readily  found  at  the  notch  from 
a  knowledge  of  the  loading  and  K   .   From  the  uniaxial  stress- 
strain  expression,  which  governs  the  ligament  at  the  tip  of  the 
notch, 

e .  £  ♦  A(ff 

We  have  shown  that 

(79) 

Putting  Equation  (78)  into  Equation  (79) ,  K   can  be  related  to 
the  local  stress  at  the  notch-tip. 

^T     Ee       EelEJ  (80) 


78) 


In  the  plastic  domain, 


.  cr 


S-Ee     K<r^ 

0  (81 
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kt  -  kj^a  (in 


Equation  (86)  can  be  written  in  terms  of  stress  concentration 
factors  to  show  the  relationship  after  yielding. 

(82) 

This  is  a  competitor  to  Neuber ' s  Relation,  which  has  been  used 
by  many  to  calculate  local  stresses  at  notches  that  are  yielded. 
A  comparison  of  the  accuracy  of  this  equation  compared  with 
Neuber' s  relation  will  be  made  later. 
RESIDUAL  STRESS 

For  each  point  along  the  loading  curve  in  the  plastic 
range,  if  unloading  occurred  at  that  point 


KT 


'.     (T-CtL 


84) 


S  (33) 

which  relates  the  modified  elastic  SCF  to  the  plastic  one. 
From  Equation  (65) , 

Kr'  =  Kr  ~  5k 
S 

This  SCF  will  be  valid  until  the  next  load  cycle  that  will  cause 
additional  yielding  is  applied. 

As  the  unloading  continues,  an  unloaded  equilibrium  point 
will  be  reached  with  a  positive  residual  strain  and  a  negative 
residual  stress.   The  value  of  the  coordinates  in  stress-strain 
space  where  this  will  occur  cannot  be  calculated  with  this  model. 
If  the  residual  strain  is  measured  photoelastically ,  the  residual 
stress  can  be  calculated  from  Equation  (65) .   The  new  stress 
concentration  factor  can  be  checked  photoelastically.   All  of 
this  has  been  done. 
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Figure  16.   Derived  Stress  -  Strain  Diagram 
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CHECK  OF  THE  THEORY 

A  70,000  lb.  load  was  applied  to  the  notched  specimen  of 
Figure  2,  which  has  a  nominal  cross-sectional  area  of  1.42  in  . 
The  nominal  stress  is  S  =  49,296  psi ,  and  the  nominal  strain 
is  e  =  4.8711  x  10~3  in/in  (E  =  10.12  x  106  psi).   By  a  pre- 
vious development  it  was  shown  that  even  for  strains  into  the 
plastic  range,  which  this  loading  does  produce  for  the  material 
at  the  notch  tip,  the  strain  concentration  factor  does  not 
change. 


e    =  KT  e  =  (2.73) (4,8711  x  10"3)  =  .0133  in/in 

Normal  incidence  photoelastic  measurements  at  the  notch 

( 8 ) 
made  by  Stenstrom     are  recorded  in  Table  VII.   Extrapolating 

to  zero  as  described  in  Part  I  gives  a  compensator  reading  of 

C   =  545.2,  or  a  fringe  value  of  11.6.   Reinforcement  effects 

of  the  coating  amount  to  a  factor  of  1.0182  when  the  material 

is  yielded,  and  the  Poisson's  Ratio  mismatch  produces  a  factor 

of  .9195.   The  actual  fringe  value  is  the  product  of  the  three, 

or  10.86.   From  the  photoelastic  fringe  equation. 


e,  -  e2  =  (1796  x  10"6)  (10.86)  =  0.01950  in/in 


At  the  notch  tip  we  are  assuming  a  uniaxial  stress  which  requires 

that  e0  =  -  ve,  where  v    =  0.5  ;  thus 
2        1 


e,  =  0.0130  in/in 
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TABLE  VII  STENSTROM  DATA 


PLASTIC  REGION  LOADING 
KT  =  2.5  Notch,  70,000  lb.  Load 


Corrected 

£1 
10"3in/in 

16.182 

£2 

10"3in/in 

-5.244 

Data  Pt. 
In. 

C 
n 

571.0 

678.0 

°1 
Psi 

163866 

°2 
psi 

0.0 

800 

0.125 

464.4 

545.8 

13.410 

-4.014 

137050 

4426 

0.250 

393.9 

472.4 

11.960 

-2.819 

125057 

12570 

0.  375 

339.7 

410.1 

10.480 

-2.264 

110367 

13360 

0.500 

294.5 

361.8 

9.476 

-1.576 

101522 

17420 

0.625 

253.3 

314.6 

8.362 

-1.143 

90502 

18160 

0.750 

219.1 

281.4 

7.807 

-0.415 

86927 

24360 

0.375 

193.0 

250.3 

7.027 

-0.215 

78828 

23710 

1.000 

170.9 

225.3 

6.436 

0.024 

73027 

24220 

RESIDUAL  COMPRESSION 


0.0 

178.0 

211.0 

-5.025 

1.654 

-50791 

43 

0.125 

158.8 

185.8 

-4.535 

1.424 

-46086 

-736 

0.250 

139.7 

163.7 

-4.004 

1.238 

-40774 

-878 

0.375 

11.7.6 

138.1 

-3.389 

1.023 

-34600 

-1016 

0.500 

100.5 

118.5 

-2.926 

0.845 

-30017 

-1316 

0.625 

81.4 

96.5 

-2.402 

0.652 

-24798 

-1551 
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This  is  a  favorable  comparison  confirming  the  uniaxial  model 
and  the  development  of  the  constancy  of  the  strain  concentra- 
tion factor.   The  error  is  2.25%. 

The  maximum  stress  is  calculated  from 

2J.SS 


Oj  ^     1,479  xio4V        3" 

1      io.iZ  *  iou  \\oaz  %  io 


Obtaining  the  value  numerically,  the  maximum  stress  is  79,403 
psi.   The  small  error  in  strain  is  insignificant  when  stresses 
are  calculated.   This  could  have  been  calculated  directly  from 
the  loading  by  using  Equation  (78) . 

The  residual  strain  was  measured  after  a  maximum  load  of 
71,000  lb.  had  been  applied.   Extrapolation  to  the  notch  tip 
gave  N_.  =  2.91,  which  determines  the  residual  strain  at 
3.48  x  10    in/in.   The  associated  residual  stress  is  obtained 
from  the  model  by  Figure  15b 


JR  "  EUR  "  £Ro: 


where  £_=£,-   4 
Ro     4    g- 

For  this  load  e.  =  0.0135  and  ?4  =  79,529  psi, 

These  values  give  eRQ   =  5.6414  x  10~3  in/in,  and  aR  =  -21,873  psi 

From  Equation  (65)  the  new  SCF  for  the  next  loading  cycle 
is  K   =2.03.   An  elastic  reloading  was  made  with  photoelastic 
measurements  taken  at  the  notch  to  evaluate  the  new  KT   .   A 
20,000  lb.  load  was  applied  which  gave  a  local  strain  reading  of 
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-3 
6.368  x  10    in/in,  which  produced  a  strain  rise  from  the  resi- 

-3 
dual  value  of  2.383  x  10    in/in  and  an  associated  stress  rise 

of  29,228  psi.   The  nominal  stress  was  20,000  lb.   acting  over 

2 

1.42  in   ,  and  the  new  SCF  was  2.08  as  compared  to  a  predicted 

value  of  2.03,  which  is  a  2.4%  error.   Several  other  specimens 
were  checked  which  resulted  in  predictions  of  similar  accuracy. 
The  actual  stress  state  for  the  above  example  is  29,228  -  21,873 
7355  psi.   This  is  a  significant  influence  on  the  fatigue  life 
calculation. 
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SUMMARY  AND  CONCLUSIONS 

A  two-dimensional  plasticity  stress-strain  law  has  been 
developed  to  coordinate  with  the  photoelastic  measurements  for 
performing  nonlinear,  plastic  analysis  of  stress  distributions 
around  notches.   These  were  compared  to  analytical  solutions 
obtained  by  slipline  theory  for  perfectly  plastic  materials. 

A  one-dimensional  model  of  the  stress  field  at  the  notch- 
tip  has  been  used  to  successfully  predict  behavior  there. 
Predictions  have  been  verified  by  experimental  measurements 
using  photoelasticity . 

The  only  remaining  piece  to  be  developed  to  have  all  the 
theory  in  place  for  fatigue  applications  is  the  one  for  cal- 
culation of  residual  stresses.   All  attempts  at  this  so  far 
have  been  unsuccessful. 
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